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Abstract 

For a quantum spin chain or ID fermionic system, we prove that the Drude weight 
D verifies the universal Luttinger liquid relation = D/k, where k is the suscepti- 
bility and Vs is the Fermi velocity. This result is proved by rigorous Renormalization 
Group methods and is true for any weakly interacting system, regardless its inte- 
grablity. This paper, combined with [T], completes the proof of the Luttinger liquid 
conjecture for such systems. 



1 Introduction and main result 

Quantum spin chains and one dimensional Fermi systems have been the subject of an 
intense theoretical investigation for decades, either for their remarkable properties or for 
the fact that they can be experimentally realized in systems like quantum spin chain 
models (KCuF3) [5] or carbon nanotubes [5]. We consider a quantum spin chain model 
with Hamiltonian 

L-1 L 

H = -J2[SlSl+, + SlSl+,]-hJ2sl + X J2 v{x-y)S',Sl + Ul, (1) 
x—i x—1 ^^^^y^L 

where S*" = c^/S for i — 1, 2, . . . , L and a 1, 2, 3, cr" being the Pauli matrices, and C/^ 
is a boundary term; finally v{x) — u(— a;) and |i'(a;)| < Ce""'^'. 

If v{x — y) = 5\x-y\^i and /i = 0, ([Ij is the hamiltonian of the XX Z spin chain in a 
zero magnetic field, which can be diagonalized by the Bethe ansatz [4 . No exact solution 
is known for more general interactions, but some particular models have been the subject 
of an extensive numerical analysis. 

It is well known that the quantum spin model can be equivalently written in terms of 
fermionic anticommuting operators = 0^=1 ("'''a)'''^! Ji = J2 = 1, one gets 

L-1 



-ff = - XI \\-'^ta^+i + 0't+i<^x]-h'^{ata^ ~ 

X—1 X—1 



2' 

l<x,y<L 
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where {/| is the boundary term in the new variables. We choose it so that the fermionic 
Hamiltonian coincides with the Hamiltonian of a fermion system on the lattice with peri- 
odic boundary conditions. 

If Ox is a local monomial in the S"" or operators, we call Ox = e^^°Oxe^^^° where 
X = (.xo,a;) and is the "imaginary time"; moreover, il A — Oxi • • • Ox„ 

< A >L.3= ^Tr[e-^H] ' (3) 

T being the time order product, denotes its expectation in the grand canonical ensemble, 
while < A >T;L,i3 denotes the corresponding truncated expectation. We will use also the 
notation < A >t= \imL,i3-¥oo < A >t;L,p- 

The response functions measure the response of the system to an external probe. In 
particular, the spin conductivity properties of model ^ can be obtained in the model ([2]) 
from the current- current response function, whose Fourier transform is defined as 



(p)= lini lim / dxo Ve*P''(Jx^o)T.L/3, (4) 

/3^oo L^oo 7_^/2 ^ ■ 



where p = {po,p), Po ^ ^n,p = ^m, (n,m) e Z , -[L/2] < m < [{L - l)/2], Jx = 
e^* Ja;e~^* and Jx is the paramagnetic part of the current 

Jx = ^[at+i'^x - a+a^+i] ■ (5) 

A crucial quantity in the study of the conductivity properties is played by the Drude 
weight, defined in the following way. Let us consider the function 

5(p) = -A-G"/^(p), (6) 

where A A^; > and 

^x = -^ [atax+i + a^+io^] (7) 

is the diamagnetic part of the current, whose mean value < A^; > is indeed independent 
of X, hence it is equal to < Ht > /L, with Ht ~ J2x ^x, the value of for /i = A = 0. 
Then the Drude weight is given by 

D= lim 5(po,0). (8) 

Po->-0 

If one assumes analytic continuation in po around pq = 0, one can compute the con- 
ductivity in the linear response approximation by the Kubo formula, see e.g. [5;, that 
is 

. . hm lim £t}^^±^ . (9) 

Therefore, a nonvanishing D indicates infinite conductivity. 

Another important quantity is the susceptibility, which can be calculated, in the 
fermionic representation, in terms of the density- density response function Gp'p(x) = 
(PxPo)t; Px = a+Gx, by the equation 

«=limG0;2(0,p), (10) 
p— >o ' 

where Ggp(p) is defined analogously to (|4]). Note that, in the fermionic representation, 
K = Kcp , where Kc is the fermionic compressibility and p is the fermionic density, see 
e.g. (2.83) of [B]. 
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The large distance behavior of the response functions is given (for couphng not too 
large) by power laws with non-universal exponents depending on all details of the Hamil- 
tonian, like the form of the potential and the value of the magnetic field. Only for the 
interaction v{x — y) — 5\x-y\.i ^ solution is known by Bethe ansatz [?], if ft, = 0; by using 
this explicit solution, k and D can be computed. However, even in that case, only a single 
exponent can be calculated 7 . 

In [SI [HI Hn] rigorous RG methods have been applied to spin chains or fermionic ID 
systems, regardless their integr ability; the outcome of such analysis is that several phys- 
ical observables, and in particular the critical exponents, can be written as convergent 
series. The exponents are interaction-dependent but nevertheless verify universal model 
independent relations; if rj is the exponent of the 2-point function (see e.g. (4) of [1]), P 
is the correlation length exponent (see e.g. (11) of [T]), X+ is the density exponent (see 
e.g. (7-9) of [T]) and X- is the Copper pair exponent (see e.g. (10) of [1 ), it has been 
proved in [UlII] that, for A small enough, 

X+ = K, X^^K-K, Y~K-^ ' 2r/-i^ + i^-i-2, (11) 

where K{X) is an analytic function such that 

;^(A) = l-A^t^^ + 0(A^), (12) 
TT smpi? 

with cospF = —h — X. Note that the exact relations ([TT|) are universal, i.e. do not depend 
on the Hamiltonian details, for instance on the form of the interaction v{x), contrary to 
the function K(X), see ([T^. 

Universal relations connect also the critical exponents with the susceptibility k; in [1] 
it was proved that 

K^—. (13) 
In this paper we prove the following Theorem for the Drude weight. 

Theorem 1.1 If X is small enough, the function I?(p) defined in (0) can be written, for 
p small hut different from Q, in the form 

Dip) = ^K P\ , +Hip), (14) 

where H{p) is a continuous function, such that \H{p)\ < C|p|'', with < ^ < 1; therefore 
the Drude weight is given by 

I? = ^ (15) 

TT 

and satisfies the identity 

vI = D/k. (16) 

The validity of the relations ([TT|). ([T ^ . ([T5|) and ([T5|) is the content of the Luttinger liquid 
conjecture formulated in [12] (see also [T3j [HI ) ; given the Drude weight and the suscep- 
tibility, one can determine exactly all the exponents and the Fermi velocity. All these 
relations are true in the Luttinger model, describing interacting fermions with a relativis- 
tic linear dispersion relation and solved by bosonization '15' ; the content of the Luttinger 
liquid conjecture is that they are true also in the model ([T|), even if the exponents are 
completely different. This is by no means obvious; the exponents, k and D are non uni- 
versal functions of the interaction, and surely depend on the dispersion relation and the 
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details of the Hamiltonian. The vahdity of the conjecture was partially checked on the 
solvable XX Z chain; Vg, k, D can be computed from the Bethe ansatz solution [T^ [TB] 
and the validity of the relation ([T5|) (following from ([TSjl . p^ ) is verified. Moreover, by 
using (fTT|) . the exponents can be exactly determined from the knowledge of k and D; note 
that the value of D found in this way agrees with the one obtained in A number of 
arguments have been proposed along the years [171 118[ 119] in order to justify the validity 
of (|lip . (jl3p and (jl5p . but they rely on unproved assumptions or approximations in non 
solvable cases. 

The present paper completes the proof of the Luttinger liquid conjecture for quantum 
spin chain or ID fermionic system with generic weak short range interaction. The proof 
relies on a number of technical results previously established in [8l[9l[T0l[TTl[T]; in particular 
the present paper extends and completes the analysis of [T], which we assume the reader 
familiar with. 

The Drude weight in quantum spin chains has been the subject in recent years of 
an intense numerical investigation [221 HSl HH with the main objective of 

detecting a possible different behavior of conductivity at finite temperatures between the 
integrable and the non integrable cases; it has been conjectured that the Drude weight 
is non vanishing also at finite temperature in the integrable cases, while it is vanishing 
in non integrable systems, but the results are still controversial. Our methods for the 
calculation of the Drude weight at zero temperature can be applied either to solvable or 
non solvable systems, and we believe that an extension of these methods would allow us 
to understand also the properties of the Drude weight at non zero temperature. 



2 Ward Identities 

We shall proceed as in App. B of [1]. Let us consider the (imaginary time) conservation 
equation: 

1^ = e''-°[H,p,]e-"''» = -idi'^J^ = -z[J,,,„ - J,_i,,J , (17) 

where we have used that [H, px] = [Ht, Px], Ht being the value of iJ for /i = A = 0. This 
equation implies some exact identities involving various correlation functions, that play 
the role in the lattice models of the usual Ward Identities (WI) of continuous relativistic 
models. They are valid at any finite /? and L, but we shall use them only in the limit 
L = /3 = oo. 

We shall call G'^(x, y) =< a^a^ > the (imaginary time) Green's function, while 

Gp^(x,y, z) Pxflya^ >t and Gj'^(x,y,z) =< J^UyU^ >t 
will be the vertex functions. By using (|17l) one gets the WI 

AG-(x,y,z)=-.a«G-(x,y,z)+ ^^^^ 
+^(a;o - zo)Sx,zG'^{y, x) - S{xo - yo)Sx,yG'^{K, z) , 

where di^"^ is the lattice derivative. In the same way a WI for the density-density correla- 
tions is derived. If we define 

Gp;p(x,y) p^Py >T, G°;j(x,y) p^Jy >T, G°'5(x,y) J^Jy >t, 
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we get 

d 



d 

^G'p;j(x,y) = G^j ji^,y) + S{xo - yo) < [P{x,xo)^J(v,xo)] > ■ 



(19) 



Noting that [p(x,xo}, P(y,xo)] = 0, while 

[P{x,xo):J(y,xo)] = -'i^x,y^(x,xo) + "^^x-l^y^y^xg), (20) 

we get, using that < Ax >=< A^; >= A, 



-zpoG°;^(p)-^(l-e-^)G",;^(p) = 0, 
-zpoG°;^^(p) - .(1 - e-^nG"j'jip) = ^(1 - e-^P)A 



(21) 



Hence, by using the definition IQ, the WI (l2Tt and the fact that G'°'j(p) = G'jp(— p), we 
get 

p2G^;2(p)-4sin2(p/2)5(p)=0. (22) 

The above equation holds quite generally for fermionic lattice systems. If ^(p) and 
D{p) were continuous in p = 0, it would imply that both k and D are vanishing. In 
the case we are considering, we will see in the next section that Gp'p(p) and D{p) are 
bounded but not continuous in p = 0, which is sufficient to prove only that: 

G°;2(po,0) = 0, 5(0,p)=0. (23) 

3 Renormalization Group anaysis 

It is well known that the correlations of the quantum spin chain can be derived by the 
following Grassmann integral, see §2.1 of l8|: 



e 



Wl,ii.m(A,J,^) ^ / p(J^j)g-V(V)+B(A,J,^)+/dx[0+^-+^-Vi-] ^ (24) 



where cospi? = —A— h — Vg = vf{1 + S), and 0^ Grassmann variables, J dx is a 
shortcut for /-^/2 ^^O' Pidi') is a Grassmann Gaussian measure in the field variables 
V'x with covariance (the free propagator) given by 

^^'^"^ ^' PL^^ -zko + {vs/vF)icospp-cosk)' 

where x(0 is ^ smooth compact support function equal to if |i| > 7 > 1 and equal to 
1 for \t\ < 1, k = (fc, fco), k • X = koxo + kx, X>l./3 = 'Dl x Vp, Vl = {k ^ 2Tm/L, n € 
Z,-[L/2] < n < [(L- l)/2]}, Vp = {k^ = 2(n + l/2)7r//3, n € Z} and 



V(^) = Ay dxdyi;(x-y)'0x'0y'0y^x + ^ j dyLtp^i^^ ~ 

-5 1 cfx[cOSpF?/'^'/'x - (V'ifeiV'x + V'x V'x+ei)/2] , 
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with ei = (0, 1), w(x — y) = S{xq — yQ)v{x — y). Moreover 

J 

^ 2 



Note that, due to the presence of the ultraviolet cut-off 7*^, the Grassmann integral has a 
finite number of degree of freedom, hence it is well defined. The constant 5m = P/^/M is 
introduced in order to take correctly into account the discontinuity of the free propagator 
(^(x) at X = 0, where it has to be defined as liva^^^t^- g(0,xo); in fact our definition 
guarantees that limA/_j.oo gj\/(x) = g{'x) for x 7^ 0, while limM^oo .9m (0, 0) = 3(0,0^). 
The density and current correlations can be written in terms of functional derivatives of 
(El 



J2 

G'p;p(x,y) = lim lim lim T^L^/g.Af (A 0, 0)L-o . 

G'p'j(x, y) = lim lim lim Wl,p,m{^^^^^)\a-o^ 

G J j(x,y) = lim lim lim ^, , ... , . Wlbm{A,0,0)\. „, 



(27) 



lim lim lim 



--^.,,,m(0,J,0)|,^„ 



In [51 [HI [TU] a multiscale integration procedure combined with Ward Identities allows us to 
write the above correlations in terms of a convergent expansion; the counterterms i/, 6 are 
chosen so that pp is the Fermi momentum and Vg is the Fermi velocity. By using Theorem 
3.12 of [S], one can easily prove that A is a finite constant. However, the bounds obtained 
from the multiscale analysis for G^'^{x,y) and Gj j(x, y) are not sufficient to prove that 
their Fourier transforms are bounded around p = (0,0). In fact, by using theorem (1.5) 
of [8], we see that their non-oscillating part behaves for large |x — y| as |x — y|"^, so that 
logarithmic divergences in the Fourier transform cannot be excluded. 

In order to compute the Fourier transform of the current-current correlation we will 
follow the same strategy used in [T] for the density-density correlation. We introduce a 
continuous model with linear dispersion relation regularized by a non local fixed interac- 
tion, together with ultraviolet 7^ and an infrared 7' momentum cut-offs. The model is 
expressed in terms of the following Grassmann integral: 

gW„(J,J,0) ^ /"p M^^g-V<"'(\/^'^+E^ = ±/'ix[Z'^)J^+^^Z(^)Xc]p,.„. 

(28) 



e 



where px.w = ^x wV'x x e A and A is a square lattice of side L, whose size is of order 
7~', say 7^'/^ < L < 7^'; P2(c?'0[''^l) is the fermionic measure with propagator 

k 

where Z and c are two parameters, to be fixed later, and xi.nO^) is a cut-off function 
depending on a small positive parameter e, nonvanishing for all k and reducing, as e — > 0, 
to a compact support function equal to 1 for 7' < |k| < 7^+^ and vanishing for |k| < 7'"^ 
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or |k| > 7^+^ (its precise definition can be found in (21) of "S^); moreover, the interaction 
is 

v(^HV') = ^ E / / - yX^V'-^^-.V'--. , (so) 

where vq (x — y) is a rotational invariant potential, of the form 

^o(x - y) = -1 ^^o(p)e'P(''-'') , (31) 
p 

with |wo(p)| < Cc^^IpI, for some constants C, /x, and vo{0) = 1. We define 



G?;^p;Jx,y,z) = lirn^ ^_ lim ^ ^ >V,.( J, J, 0)|,,j,,,o , 

d ~ 
G'?AU^(x,y,z) = lim lim^ T^Tr-+-^r7— W;,w(^. j=J=0=o ' 

GLjy,z)= lim hm W,Ar(J,J,0)|,^j,^^o, (32) 

G°/',/.j(x,y) - lim lim ~ ~ >Vi,jv(J, J,0)|j^j^^^o- 

-l,N-^oo a~^,L^oo gj^QJy 

The existence of the N ^ oo limit has been proved in (26j and in §3 of [TT] , extending the 
method used in |27] for the analysis of the Yukawa model in two dimensions; the existence 
of the limit I ^ —oo has been proved in [8l l9l [TO]. 

The model is a sort of effective model for the lattice fermionic model (2); it is 
indeed well known that a non relativistic gas of fermions in one dimension admits an 
effective description in terms of massless Dirac fermions in d = 1 + 1 dimension. We can 
make precise this idea via the following lemma, whose proof is an immediate extension of 
the proof given in §3 of [IJ for the density-density correlation. 

Lemma 3.1 Given X small enough, there exist constants Z, Z^^\ Z^^\ Aqo, depending 
analytically on X, such that Z = 1 + ©(A^), Z^^) = l + 0(A), Z^^) = vp + 0(A), A^o = 
A + O(A^) and, if c = Vg and |p| < k < 1, 

with Ap^p{p), Aj^jij)) Lipschitz continuous in p. Moreover, if we put p'^ — (0,u;pi?) and 
we suppose that < k < |p|, |k'|, |k' — p| < 2k, < d < 1, then 

Gl'^ik' + , k' + p + p^) = GfilJk', k' + p)[l + 0{k^)] , 
GY{k' + p-p,-k'+p + p'^p) = Gflj.Jk',-k' + p)[l + 0{K^)] , (34) 
G\k' + p-p) = GlJk')[l + OiK^)]. 

This lemma says that the vertex functions of the two models are essentially coinciding 
close to the Fermi momenta, if the bare parameters are chosen properly, while the response 
functions differ by a continuous function. Note also the the bare parameters of the model 



d 
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are expressed by convergent expansions depending on all model details, but the WI 
imply that they are not independent parameters, as we will see shortly. 

The main reason behind the introduction of the model ([25)1 is that, while the model 
ID) is invariant only under the phase transformation tp^ e^'^"'ijj^, the model ([25)) is 
invariant under two phase transformations, the total "^xw ^ c±«a,/,± 



e^""V'x w the chiral 

V'x w This implies that the Fourier transforms of the response functions can 
be completely determined from the WI, see app. A of 1 ; if D^{p) = —ipo + ojcp, we get: 



■0,2 
th,J,J 



0,2 



-1 

47rcZ2 1 - r2 

-1 



D^ip) , D+{p) 



^+(P) 



^-(P) 
^+(P) 



D+ip) D^{p) 



2r 



2r 



0(p) 
0(p) 



(35) 



where r — j^. Therefore, from ([55]) and (j33p . since c ^ 



G°;^(P) = 

D{p) = 



-1 



2 r 



D^(p) D+jp) 
D+{p) D^ip) 

D^ip) , D+{p) 



2t 



Ap,p{0) + Rp{p), 



D+ip) D^(p) 



- 2t 



Rj{p), 



A 



(36) 



with \Rp{p)\,\Rj{p)\ < C\p\^, < 1? < 1. The constants Ap,p(0), Ajj{0) and A are 
expressed by convergent expansions, but their values can be determined from the WI for 
the model ([1]); indeed, by ([55)) . G'^j^j j and G^/f. „ are not continuous in p = 0, but they 



are bounded, so that ([25)) holds; this condition fixes the values of Aj.j(O) - 
so that 



-A and Ap^p{0) 



1 



Z(3))2 



9 2 



5^2 1 - r2 p2 _^ z;2p2 



1 (Z(3))2 



r2 p2 



y2p2 



hi?j(p). 



(37) 



Moreover, the vertex functions verify the following WI, see (35) of [1]: 



G?AL^(k,k + p)+a;pi;, 



1 



-G 



1 



[G2 (k)-G2 



th-A 



(k + p)]; 



2,1 



.(k,k + p) = 



(38) 



hence, by using (l34l) and by comparing ([38)) with the WI 
parameters are not independent, but verify the relations: 



implying that 



(l-r)Z 



= 1, 



^(3) 

Vs — = 1 , 

Z(3) 



^PP{P) = 



K 



9 7 
VsP 



D{p) = ^K— 

7^ Po 



h f 2p2 



+ Rp{p): 



^ +i?j(p). 



18]), we get that the bare 



(39) 



(40) 



with K — jzfT- Eq. (52) of 8 shows that K is indeed the critical index see (ITT)) : 
hence, by using ^ and (flUl) , we get the relations ([T5)) and (IT51) , which immediately imply 
(Unj), so that Theorem 11.11 is proved. 
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